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ABSTRACT 


This  paper  is  an  analysis 
shot  KH261  when  fired  from  th 
Density  shot  ammunition  is  a 


of  the  behavior  of  density 
e  .45  caliber  piste  1  {M1911A 
multi-projected  round  which 
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contains  16  snail  pellets.  Both  the  standard  rifled  barrel 
s  ijtsiCOti  o  a r  r  e  1  ^ere  us  wi  » o  to  llect  data  on  the  behavior 
of  the  XM261  round.  A  family  of  distributions  is  selected 
for  each  barrel  type  to  model  the  scatter  of  pellets  from 
the  XM2S1  round  at  varying  ranges. 
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III.  SELECTION  OF  DISTRIBUTIONS 


For  the  present  analysis,  the  basis  for  selecting 
distributions  to  model  the  scatter  of  pellet  strikes  from 
the  XM261  was  a  letter  from  the  United  States  Army  Material 
Systems  Analysis  Agency  (USAMSAA)  [Ref.  2],  USAMSAA  was 
replying  to  a  request  from  USASASA  to  analyze  data  obtained 
from  firing  the  XM261  round.  The  data  analyzed  by  USAMSAA 
was  prepared  by  the  Frankford  Arsenal.  As  a  result  of  their 
analysis,  USAMSAA  suggested  that  the  distribution  of  pellet 
strikes  from  the  XM261  round  when  fired  from  a  rifled  barrel 
could  be  modeled  by  using  a  polar  coordinate  system  with  the 
angle  having  a  uniform  distribution,  the  radius  vector  having 
a  normal  distribution,  and  the  angle  and  radius  vector  being 
independent.  (The  model  based  on  these  assumptions  will 
hereafter  be  called  the  Polar  Model.)  When  the  XM261  round 
was  fired  from  a  smooth  bore  barrel,  USAMSAA  suggested  that  a 
bivariate  normal  distribution  with  independent  X  and  Y  would 
serve  as  an  adequate  model.  This  student  felt  that  these  two 
models  could  serve  as  an  adequate  basis  for  an  initial  investi¬ 
gation  of  the  XM261  data  provided  by  H.  P.  White  Laboratory. 

The  plane  had  to  be  partitioned  into  regions  so  that  a 

chi-square  goodness-of-f it  test  could  be  performed.  The 

objective  for  the  partition  was  to  allow  the  maximum  number  of 
« 

regions  consistent  with  the  requirement  to  have  an  expected 
frequency  of  observation  of  pellets  of  approximately  five 
[Ref.  3  and  4],  In  connection  with  testing  the  fit  of  the 


10 


Polax’  Modal ,  the  plane  was  originally  divided  into  30  nearly 
equally  probable  regions.  Since  the  angle  was  assumed  to 
have  a  uniform  distribution,  the  plane  was  first  divided  into 
six  regions  by  constructing  six  rays  from  the  origin  at  50 
degree  intervals.  Then  four  circles  were  superimposed  over 
the  rays  to  further  divide  the  plane.  The  radii  of  the 
circles  corresponded  to  the  20,  40,  60,  and  80  percentiles 
of  the  fitted  distribution  of  the  radius  vector.  This 
resulted  in  a  parti tion  of  the  plane  into  30  nearly  equally 
probable  regions  (see  Figure  1) ,  Since  there  were  160  data 
points,  the  expected  frequency  for  each  of  ths  regions  was 

5.33. 

When  testing  the  fit  of  the  b  -~.riate  normal  distribution, 
the  plane  was  divided  into  36  rectangles.  The  mean  and 
standard  deviation  of  the  X  variable  were  estimated  using 
the  maximum  likelihood  estimates.  Then  values  for  the  16.67, 

33.33,  50,  66.67,  and  83.33  percentiles  of  the  X  variable 
were  estimated  using  the  estimates  for  the  mean  and  standard 
deviation.  A  vertical  line  was  constructed  on  the  X  value 
which  corresponded  to  each  ot  these  estimated  percentiles. 
Then  using  the  same  method  horizontal  lines  were  constructed 
on  the  same  sample  percentile^  of  the  Y  variable.  The 
resulting  partition  provided  36  nearly  equally  probable 
rectangles  (see  Figure  2) .  The  expected  frequency  based  on 
the  160  data  points  was  4.45.  Although  this  was  less  than 
the  desired  frequency  of  five,  it  was  sufficient  [Ref.3j. 
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IV,  ANALYSIS  OF  DATA  FROM  THE  RIFLED  BARREL 

A.  INITIAL  RESULTS 

When  fitting  tha  bivariate  normal  distribution,  four 
parameters,  the  mean  and  standard  deviation  of  both  the 
X  ana  Y  variables,  had  to  be  estimated.  Since  the  X  and  Y 
variables  were  assumed  to  be  independent,  there  was  no  need 
to  estimate  a  correlation  coefficient.  Because  there  were 
36  regions  and  four  parameters  estimated,  this  chi-square 
test  had  31  degrees  of  freedom.  This  yielded  a  critical 
value  of  45.0  for  a  fi/e  percent  level  of  significance  and 
a  critical  value  of  41.5  for  a  ten  percent  level  of  signifi¬ 
cance.  The  test  statistic  had  a  value  of  95.6  for  the  10 
meter  data  and  a  value  of  93.8  for  the  30  meter  data. 

Hence,  the  hypothesis  that  either  set  of  data  was  distributed 
according  to  a  bivarate  normal  was  rejected. 

When  fitting  the  Polar  Model,  only  two  parameters,  the 
mean  and  standard  deviation  of  the  distribution  of  the 
radius  vector,  were  estimated.  Eecause  the  Polar  Model  had 
30  regions  and  two  parameters  were  estimated,  this  chi-square 
statistic  had  27  degrees  of  freedom.  The  critical  values  for 
five  and  ten  percent  levels  of  significance  are  40,1  and 
36.7,  respectively.  The  test  statistic  for  the  10  meter  data 
was  55.25,  that  for  the  30  meter  data  was  also  55.25.  This 
led  to  the  rejection  of  the  hypothesis  that  either  set  of 
data  was  distributed  according  to  the  Polar  Model . 
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Since  both  distributions  had  been  rejected  as  methods 
for  adequately  describing  the  data,  an  additional  distribution 
war  sought,  A  histogram  of  the  values  for  the  radius  vector 
suggested  that  the  radius  vector  might  be  distributed 
according  to  a  gamma  distribution.  Therefore#  an  additional 
distribution  was  introduced.  This  distribution  utilised  the 
polar  coordinate  description  of  the  data  with  the  assumptions 
that  the  angle  was  uniformly  distributed  between  zero  and 
two  pi#  the  radius  vector  was  distributed  according  to  a 
gamma  distribution#  and  the  angle  and  the  radius  vector 


were  independent.  This  distribution  will  be  referred  to  as 
the  Polar  Gamma  distribution.  The  plane  was  partitioned  for 
the  chi-square  goodness-of-f it  test  in  the  same  way  as  for 
the  Polar  Model.  The  values  of  the  radii  of  the  circles 
corresponded  to  the  appropriate  percentiles  of  the  gamma 
distribution  of  the  radius  vector. 


If  a  gamma  distribution  was  to  be  used  to  represent  the 
distribution  of  the  radius  vector#  the  parameters  would 
have  to  be  estimated  from  the  data.  The  maximum  liklihood 
method  of  estimating  the  parameters  has  some  appeal,  but 
it  is  cumbersome  because  it  involves  the  gamma  function. 
Therefore,  it  was  decided  to  use  the  method  of  moments  to 
estimate  the  gamma  parameters..  This  decision  was  made  in 
spite  of  the  fact  then,  me  method  of  moments  produced  a 
biased  estimate  of  the  parameters.  Appendix  B  contains  the 


method  of  moments  estimates  and  a  justification  that  they 
are  biased.  Even  after  the  parameters  were  estimated  some 


difficult f  is  experienced  in  estimating  the  20,  40,  50, 
and  80  percentiles  of  a  general  gamma  distribution  to  use  a 
values  for  the  radii  of  the  circles  separating  the  different 
regions  it  the  plane.  However,  if  the  shape  parameter 
alpha  is  a  positive  integer,  then  a  convenient  relationship 
exists  between  the  gamma  and  the  Poisson  distributions. 
Specifically,  if  X  is  distributed  gamma  {M, J)  and  Y  is 
distributed  Poisson  (Xx)  t  then  P  {X£%)=P  C'Y>K)  [Ref.  5] . 
Therefore,  the  percentile  values  for  the  gamma  distribution 
with  an  integer  alpha  value  can  be  interpolated  from  a  table 
of  cumulative  values  for  the  Poisson  distribution. 

One  of  the  considerations  for  this  problem  was  that 
the  distribution  used  to  describe  the  data  should  lend 
itself  to  easy  generation  by  computer  simulation.  Hence , 
it  was  decided  that  after  the  gamma  parameters  were  estimated 
using  the  method  of  moments,  the  integers  on  either  side  of 
the  estimated  alpha  value  would  be  usee  instead  of  the 
estimated  alpha  value.  This  approach  would  permit  the  use 
of  the  relationship  between  the  gamma  and  the  Poisson 
distributions ,  Each  data  set  would  then  yield  tvo  chi-square 
test  statistics  for  each  set  of  gamma  parameters  estimated; 
One  test  would  use  the  integer  value  immediately  below  the 
estimated  alpha  value,  the  second  would  use  the  integer 
value  immediately  above  the  estimated  alpha  value.  Both 
tests  would  use  the  estimated  lambda  value. 


Since  fitting  the  polar  gamma  also  requires  the  estimatic 
of  two  parameters  and  the  piano  had  been  partitioned  into  30 


regions ,  tne  chi-square  goodness -or "fit  test  ene  Poxar 
Gamma  had  27  degrees  of  freedom.  Consequently*  the  critical 
values  for  the  test  statistic  were  the  sane  as  for  the  Polar 
Model,  The  values  for  the  two  test  statistics  generated 
using  the  10  meter  data  with  the  rifled  barrel  were  62,0  and 
4a,l.  Those  for  the  30  meter  data  with  the  rifled  barrel 
were  58.6  and  52.6*  Therefore,  both  of  the  null  hypotheses 
that  the  two  sets  of  data  were  distributed  according  to  the 
Polar  Gamma  distribution  were  rejected.  Although  the 
bivariate  normal.  Polar,  and  Polar  Gamma  distributions  had 
boen  rejected  as  models  for  the  data,  the  Polar  Gamma 
distribution  had  produced  the  lowest  value  for  the  test 
statistic  for  both  the  10  and  30  meter  data.  This  indicated 
that  the  Polar  Gamma  distribution  might  have  some  potential 
to  describe  the  data. 

Because  all  the  distributions  had  failed  to  produce 
an  acceptable  description  of  the  data,  a  check  was  made 
to  see  if  the  data  contained  any  identifiable  variation 
that  had  not  been  accounted  for.  The  physical  character¬ 
istics  of  the  method  of  data  collection  led  to  the  belief 
that  such  a  variation  might  exist.  The  recoil  and  rifling 
of  the  barrel  were  suspected  of  causing  +  he  center  of  mass 
{sample  mean)  of  a  pellet  pattern  to  drift  up  and  to  the 
left  as  the  range  increased.  The  recoil  was  suspected  of 
causing  the  drift  up.  The  movement  -i  the  barrel  of  one 
half  of  one  degree  is  sufficient  to  cause  the  center  of  mass 
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inches  at  30  meters .  The 


fjung  was 
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suspected  o;  causing  the  drift  to  the  left.  Since  the 

left  handed  rifling  produces  a  counterclockwise  rotation 

of  the  projectile,  the  project  Me  would  tend  to  draft  to  the 

left  when  meeting  air  resistance.  If  the  center  of  mass  did 

in  fact  drift  from  the  aim  point  and  if  the  drift  was  random, 

then  this  variation  in  the  location  or  the  centers  or  mass 

would  introduce  an  uncontrolled  variation  among  the  rounds. 

This  would  be  considered  undesirable  since  the  ehi-sguare 

goodness-o£-£it  test  assumes  that  all  observations  tshcts) 

come  from  the  same  distribution. 

In  order  to  determine  whether  the  data  would  tend  to 

support  or  refute  the  hypothesis  that  the  centers  of  mass 

did  drift  as  range  increased,  the  center  of  mass  of  each 

shot  was  calculated  from  the  data  at  ranges  of  10  and  30 

meters .  Each  center  of  mass  was  based  cm  the  16  pellets 

in  one  Xii2ci  rr.  and.  hn  examination  of  the  centers  of  mass 

calculated  in  this  fashion  tended  to  support  the  hypothesis 

that  the  center  of  mass  drifted  up  and  to  the  left  as 

range  increased;  For  both  the  10  and  30  merer  data  eight 

out  of  ten  X  coordinates  of  the  centers  of  mass  were  negative 

and  nine  out  of  ten  Y  coordinates  were  positive.  Furthermore , 

a  comparison  of  the  centers  of  mass  at  10  and  30  meters 

for  the  same  round  revealed  that  eight  out  of  ten  centers 

of  mass  had  drifted  up  and  to  the  left.  It  was  also  noted 

that  the  actual  distance  that  the  center  of  mass  moved 
■% 

varied  significantly  among  the  rounds.  The  calculated  cen¬ 
ters  of  mass  are  provided  in  table  1. 
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I 


Overall  Mean 


TABLE  I 


CENTERS  OF  MASS-RIPL-E.0  BARREL 
10  Meters 
(0.25,1.14) 

C“ 0.10,0-75? 
(-0.41,2.88) 

(0 .00 ,1.35) 
(-1.92,-1.31) 

(-0 .43,2.33) 

(-0.85,0.64) 

(-0.53,1.48) 

(-2.50,2.27) 

(1.29,0.05) 

(-0.52,1.16) 


30  Meters 
(-0.46,2.25) 
(-1.27,2.30) 
(0.50,8.35) 
(-3. 03, 2. S3) 
(-3,76,-4.30) 
(-8.43,3.92) 
(-1,45,2.15) 
(-0.88,2.51) 
(-5,37,2.58) 
C-4.08,0.56) 
(-2.62,2.36) 


the  centers 


center  or  aass  at  tne  origin. 


An  interesting  phenomenon  may 


:ur  wnsn  eacn 


10  of 


pe-J Vats  has  a  center  of  mass  at  the  origin, 


The  data  raaj 


appear  to  have  an  area  of  relatively  low  probability  at  the 
origin  {i.e«,  a  ’hole'  may  appear  in  che  data) .  This  may  be 
true  even  though  the  original  data  may  not  have  exhibited 
this  trait.  A  possible  explanation  for  the  ‘hole*  appearing 
is  that  when  the  centers  of  mass  are  not  constrained  to  be 
at  the  origin,  the  area  of  lew  probability  :or  one  shot  may 
lie  in  an  area  of  high  probability  for  a  different  shot. 

The  overlap  of  the  different  patterns  from  trae  different 
rourfe  may  prevent  a  ‘hole’  from  being  observed.  However , 
when  the  centers  of  mass  are  all  at  the  origin,  the  areas 
of  lew  probability  may  all  coincide  to  produce  a  'hole*  in 
the  data. 

Another  phenomenon  may  produce  a  ’hole'  in  the  data  even 
when  the  distribution  of  the  data  does  not  exhibit  an  area 
of  low  probability.  For  example,  if  one  drew  pairs  of 
numbers  from  a  uniform  distribution  and  plotted  the  distance 
of  each  paint  from  the  center  of  mss  of  its  pair  of  points, 
then  he  would  observe  an  area  of  low  probability  at  the 


r 


Although  the  modification  of  the  data  seemed  convenient 
to  allow  for  its  analysis ,  this  procedure  introduced  three 
significant  problems .  First ,  it  complicated  the  problem  of 


utilizing  any  distribution  to  measure  the  effectiveness  of 
the  round  by  US AS AS A.  Any  distribution  produced  by  subs  ?quent 
analysis  would  reflect  the  scatter  of  pellets  about  the 
center  of  mass,  but  would  suppress  the  drift  of  the  center 
of  mass  from  the  aim  point.  The  drift  of  the  center  of  mass 
would  have  to  be  reintroduced  in  ar:y  simulation  of  the 
actions  of  this  round.  Failure  to  do  so  might  produce 
estimates  of  the  effectiveness  of  this  round. 


erroneous 


ths  distribution  or  toe  eata  ...uen  a  Divariate  ootaaj. 
distribution  is  used#  one  would  suspect  that  a  similar 
phenomenon  might  occur  with  other  distributions.  Hence, 
centering  the  pellets  at  their  respective  center  of  mass 
say  appear  to  increase  the  probability  of  striking  close 
to  the  origin. 

The  final  problem  produced  by  modifying  the  data  concern: 
the  degree  to  which  the  data  conforms  to  the  assumptions 
required  for  the  chi-square  coodness-of-fit  test  after  xt 
has  been  modified.  Even  if  the  data  bad  been  independent 
in  its  original  fora,  it  is  important  to  note  that  after  the 
data  has  been  modified,  the  coordinates  of  the  pellets  in  on 


two  tolar  Gazr.a  Models  could  net  be  re^ected- 

A  similar  pattern  emerged  - her.  toe  30  sister  data  was 
analyzed.  Use  test  statistic  with  tbs  bivariate  noraal 
distribution  was  83.4.  That  for  the  Polar  Model  was  30. S 
the  Polar  Gamma  distribution  with  alpha  equal  to  15  and 
iamda  equal  to  0.497/  yielded  a  test  statistic  of  30.5. 
Wien  the  gamma  parameters  were  alpha  equal  to  16  and  lams 
equal  to  0.4377,  a  test  statistic  of  27.1  was  obtained. 


Aerain  the  Polar  Jfode 
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two 
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distribution.  For  the  30  meter  data  the  Polar  Gamma  distri¬ 
bution  with  alpha  parameter  equal  to  16  would  he  preferred 
above  all  other  distributions.  However,  when  choosing 
between  the  two  Polar  Gamma  distribution  for  the  10  meter 
data,  an  additional  consideration  should  be  included.  It 
seemed  reasonable  that  the  distribution  selected  to  model 
the  scatter  of  pellets  at  10  meters  should  resemble  the  30 
meter  distribution  as  closely  as  possible.  The  addition  of 
the  consideration  favors  the  selection  of  the  Polar  Gamma 
distribution  with  alpha  equal  to  16,  even  though  the  test 
statistic  for  the  other  Polar  Gamma  distribution  is  slightly 
smaller.  Therefore ,  this  student  would  select  the  Polar 
Gamma  distribution  with  alpha  equal  to  16  as  a  model  for  the 
scatter  of  pellets  fired  from  a  rifled  barrel.  Additional 
justification  for  preferring  the  Polar  Gamma  distribution 
over  the  Polar  distribution  is  found  in  the  characteristics 
of  the  different  distribution.  Since  the  Gamma  distribution 
is  defined  only  for  positive  values,  using  a  Gamma  distri¬ 
bution  to  mode}  the  distribution  of  values  of  the  radius 
vector  is  appealing.  The  normal  distribution,  on  the  other 
hand,  is  unbounded  in  both  the  positive  and  negative  direc¬ 
tions.  -This  always  allows  the  theoretical  possibility  of 
having  a  negative  value  for  the  radius  vector.  If  the  mean 
of  the  distribution  of  the  radius  vector  is  not  large  with 
respect  to  the  standard  deviation,  using  a  normal  di;.cribu- 
tion  could  produce  a  significant  number  of  negative  values. 


Moreover/  the  histogram  of  data  values  for  the  radius 
vector  appeared  to  resemble  a  gamma  distribution  more  than 
a  normal  distribution.  Finally,  the  gamma  distribution  is 
easily  modified  to  allow  for  variation  in  ranges.  Only 
the  scale  parameter  lambda  need  be  modified  to  reflect  a 
change  in  range.  Possible  procedures  to  account  for  the 
change  in  lambda  as  the  result  of  a  change  in  range  are 
discussed  in  appendix  F. 


V.  ANALYSIS  OF  SMOOTH  BARREL  DATA 


The  analysis  of  the  smooth  barrel  data  was  similar  to 
the  analysis  of  the  rifled  barrel  data.  The  parameters  for 
thfe  Gamma  distribution  were  estimated  using  the  method  of 
Elements.  The  three  distributions,  bivariate  normal,  Polar, 
and  Polar  Gamma,  were  then  used  to  determine  if  they  could 
adequately  explain  the  data.  All  three  failed  to  provide  an 
acceptably  small  value  for  the  test  statistic.  For  example, 
the  10  meter  data  produced  test  statistics  of  145.5,  291.0, 
359*2,  and  248,5  for  the  bivariate  normal,  polar,  and  the 
two  versions  of  the  Polar  Gamma  distributions,  respectively. 

Following  the  technique  of  the  previous  section,  the 
centers  of  mass  were  calculated  using  the  data  points  in 
each  round.  These  calculated  centers  of  mass  were  then 
examined.  As  with  the  rifled  barrel,  the  smooth  bore  data 
centers  exhibited  a  random  drift.  However,  the  characteris¬ 
tics  of  the  drift  ware  slightly  different.  Since  the  barrel 
was  no  longer  rifled  and  since  the  rifling  was  suspected  of 
producing  the  drift  toward  the  left,  one  would  expect  the 
smooth  barrel  to  fail  to  consistently  produce  a  movement  to 
the  left.  The  smooth  barrel  data  was,  in  fact,  free  from 
this  tendency.  Essentially  half  cf  the  centers  of  mass 
drifted  to  the  left  and  half  to  the  right.  There  continued 
to  be  a  tendency  for  the  centers  of  mass  to  rise  as  the 
range  increased.  As  a  matter  of  fact,  all  ten  centers  of. 
mass  rose  as  the  range  increased.  In  order  to  remove  this 


random  drift,  the  data  was  again  modified  by  subtracting  the 
coordinates  of  tie  center  of  mass  from  the  coordinates  cf 
each  impact  poin in  that  round.  The  centers  of  mass  for 
the  smooth  barrel  data  are  provided  in  Table  2. 

However,  evon  with  the  data  modified  in  this  way ,  the 
values  of  the  last  statistic  were  unacceptably  large.  Only 
the  Polar  Gamma  distribution  produced  a  test  statistic  less 
than  100.  While  this  student  was  preparing  a  histogram  of 
values  for  the  radius  vector,  he  noticed  that  rounds  6,  7, 
and  8  produce!  an  inordinate  number  of  large  values  for  the 
radiua  vector.  Ibis  caused  hist-  to  suspect,  that  these  rounds 
might  be  classified  as  outliers,  Consequently,  he  caicula** 
ted  the  sample  variance  for  the  three  rounds  in  question  and 
the  sample  variance  for  the  remaining  seven  rounds.  When 
the  ratio  of  these  sample  variances  was  calculated,  an  F 
statistic  was  produced.  The  F  statistic  for  tne  10  meter 
data  was  21.9,  for  the  30  meter  data  30.3.  The  critical 
value  f t r  a  level  of  significance  of  five  percent  and 
degrees  of  freedom  46,112  is  approximately  1.55.  Since  the 
test  statistics  were  large  with  respect  to  the  critical 
value,  this  student  classified  rounds  6,  7,  and  P  outliers 
and  removed  them  from  the  data  set.  This  variation  among 
the  rounds  may  have  been  produced  by  variation  in  the  manu¬ 
facture  or  performance  of  the  ammunition.  In  any  event, 
there  setras  to  be  sufficient  reason  to  desire  a  repeat  of 
this  experiment  to  see  if  this  variation  is  observed  again. 


If  30  regions  had  been  used  for  the  goodness -of -fit  test 
after  the  data  set  had  been  reduced  by  the  exclusion  of  the 
three  outlying  rounds,  this  would  have  resulted  in  an  expected 
frequency  of  only  3.73,  much  less  than  the  desired  five. 


Consequently,  the  number  of  regions  used  when  testing 


the  fit 


of  the  Polar  and  Polar  Gamma  distributions  was  reduced  to  24 
by  superimposing  only  three  circles  over  the  six  rave  from  the 
origin.  The  radii  corresponded  to  the  estimated  25,  50,  and 
75  percentiles  of  the  distribution  of  the  radius  vector.  The 
use  of  24  regions  raised  the  expected  frequency  to  4.67 
which  is  acceptable  fPef.  3} . 

After  eliminating  the  atypical  rounds,  reducing  the  number 
of  regions  in  the  partition  of  the  plane,  and  reestisating 
the  gaicsa  parameters  based  on  the  remaining  seven  rounds, 
the  three  dis  ributions  were  tested  for  fit.  The  bivariate 
normal  a&d  Polar  distributions  resulted  in  relatively  large 
test  statistic sj  all  were  over  30.  The  Polar  Gamma  produced 
a  test  statistic  PF  23.3  for  parameter  values  of  alpha 
equal  to  4  and  i arnica  equal  to  1.8934,  using  the  10  meter  data. 
The  lowest  test  statistic  for  the  30  meter  data  was  24.7, 
obtained  with  the  Polar  Gamma  distribution  with  parameter 
estimates  of  alpha  eoua_  to  4  and  lambda  equal  to  0.5352. 

Since  u.v.  smooth  barrel  data  had  been  modified  ir,  the  same 
manner  as  the  rifled  barrel  data,  the  same  problems  were  encoun 
iered.  Since  the  Polar  Gamma  distribution  with  an  alpha  value 
of  4  produced  the  lowest  value  of  the  test  statistic  for 
both  sets  of  data,  it  was  recommenced  to  model  the  behavior  of 


the  smooth  barrel  data 


>.s  was  tne  case  with  the  rifled  data,  the 
distribution  was  preferred  over  the  Tolar  dis 
of  the  ehacteristics  of  the  gasaa  distributee 
histogram  of  values  for  the  radius  vector  spo 
gam&a,  and  because  the  Polar  Gamma  lent  icsel 


Polar  Gamma 
fribution  because 
n ,  because  the 
eared  to  be 
£  to  modification 


as  the  ranoe  varied 


Shot 


APPENDIX  A 
ACCUFJsCX  DAvA 
Smooth  Bore  Barrel 


Impact 

at  10  Meters- 

Xisp3Ct  d  ts 

30  aet< 

XCin.) 

'l  (in.) 

X{in.) 

YCin 

-0,1 

+0. 4 

+3,2 

-2.5 

4-0,4 

+0.5 

-2.4 

—2.6 

+0.6 

+3.3 

-1.8 

-1,7 

+„.S 

— 3 .  < 

+0.1 

—0 ,8 

+1.0 

+5.8 

+2.2 

-0.4 

+1.3 

+2.0 

+2 . 4 

+1.1 

+K4 

+2,8 

-0.3 

+2.3 

+0.4 

+3.7 

-1.0 

+2,? 

-2.8 

+4  .7 

-0.8 

+2.8 

-1.4 

+4,8 

+2.4 

s 

+8.1 

+7,4 

+2,0 

+3.3 

-2.1 

*7,5 

+1.7 

+3.4 

+  -0  9 

+7.7 

-3,9 

+3!? 

■«+:  4 

-v  *  - 

t/.§ 

-0.1 

+3.8 

AA  +  W 

+S.I 

+0.4 

+4.1 

+2.6 

+9.1 

+1.1 

-0.1 

+3.4 

-3.7 

0.4 

*0  .  9 

-1-2 

-0.5 

-0.1 

+1.3 

-3,1 

-0.5 

-0.6 

+1.8 

0.0 

+0.3 

+0.9 

+2.6 

+2-,^ 

+2.9 

-1.4 

+2.7 

-4.5 

+3.9 

+0.S 

A/  *| 

5  *  *  * 

+0.9 

+5.4 

+0.9 

i’S  ^ 

+1.4 

+5.5 

+4.1 

+3.9 

-4  2 
—  « 

+7.4 

+0 . 5 

+4.0 

^1  1  £ 

—  1  ,v 

+0.6 

+4,4 

+0.3 

+8.0 

-0,7 

•  *  * 

— -s  .  ** 

■f§.4 

+8.7 

-0,8 

+4.7 

—  *  ” 

+10.2 

+0.8 

+5,6 

+5.9 

+10.6 

+0.9 

j-5.9 

+1.1 
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-2.1 

+€  .7 
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+16 . 6 
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Shot 

Impact 

at  10  Meters 

Impact 

at  30  Meters 

Nurrbea; 

X(in.) 

Y (in.) 

X(in. ) 

Y  (in. 

3a 

+3,4 

-7.4 

+35.3 

+6.3 

b 

-0.9 

-8.6 

+  33.4 

+  33.5 

c 

-1.8 

-8.1 

+33.0 

+33.4 

d 

-3.7 

-S.5 

+18 .  j 

+36.5 

e 

-7.4 

+1.8 

+8.1 

+38.4 

f 

-6.5 

+3.2 

+10.3 

+4.8 

g 

-11.3 

+1.8 

+2.7 

+34.3 

h 

-12.6 

+3.7 

-1.0 

+40.0 

i 

-5,2 

+5 . 4 

-21.4 

+5.5 

j 

-3.6 

+13.1 

-23.1 

+2.5 

k 

+0.9 

+13.8 

-37.2 

+1 . 3 

1 

+1 . 8 

+15.7 

-40.2 

+7.2 

m 

+2.2 

+15.7 

-2.7 

-28.5 

a 

+4.0 

+1.9 

-4.9 

-27.0 

o 

+12.3 

+3.0 

-11.6 

-31-6 

P 

+14.6 

+4.3 

+9.1 

-23,1 

4a 

412.4 

-7.G 

+29.5 

+32.4 

b 

+12.2 

-6 , 6 

+24.0 

+35.7 

c 

-;7.o 

-10.0 

+21.5 

+29.5 

d 

+7.5 

-4.6 

+16 . 2 

+15.9 

e 

-4.0 

-9.4 

+38.4 

-19.0 

f 

-6.4 

-5.9 

+35.1 

-16.0 

<3 

-8.0 

-6.6 

+21.0 

-12.3 

h 

-10.9 

-7,5 

+6 . 1 

-10,1 

i 

-9,0 

+7.4 

-7.5 

-32.3 

j 

-9.9 

+9.8 

-25.2 

-23.8 

k 

-4.7 

+12.2 

-38.8 

-26.5 

1 

-13.3 

+14.4 

-30.3 

-39.7 

m 

+8.4 

+11.1 

-31.2 

+21.6 

n 

** 

+8.0' 

+9.5 

-33.5 

+28.5 

o 

+10 . 3 

+10.3 

-33.5 

+  26,2 

p 

T#  ♦  i 

+4  ,7 

-38.4 

+35.1 

Sequence  of  sub-projectiles  lister!  at  2  -  niters  is  not 
necessarily  identical  with  the  secruence  listed  at  10  meters 


Shot 

Impact 

at  10  Meters 

Impact 

at  30  Meters 

Number 

X(in„) 

x (in.) 

X{in. ) 

V (in. ) 

5a 

+12.8 

-13.6 

+12.0 

+  17.2 

b 

+7. 

-4.1 

+2.5 

+24.0 

c 

+1.4 

-2.5 

+18.6 

+0.1 

d 

+7.3 

+0.2 

+14.0 

+21.5 

e 

-3.0 

-11.5 

+2,5 

-6.7 

f 

-11.7 

-15.6 

+22.5 

-9.5 

g 

-4.8 

-8.3 

+40.2 

-39.5 

h 

-15.3 

-5.9 

+13.8 

-35.5 

i 

-14.8 

-2.0 

-4.4 

-33.5 

3 

-15.7 

-2.2 

-14,2 

-37.0 

k 

-12.3 

+1.5 

-14.1 

-24.0 

i 

-0.1 

+13.0 

-35.8 

-47.5 

m 

+4.2 

+15.3 

-42.0 

+7,7 

n 

+3.6 

+7.3 

-30.5 

+  32,5 

o 

+6,0 

+5.4 

-26.7 

+30.2 

P 

+6.0 

+2.Q 

—23 . 6 

+31.2 

6a 

+11.5 

+1 . 5 

+0.6 

+34.6 

b 

+9.5 

-0.9 

*2.4 

+37.6 

c 

+7.9 

+0.3 

+23.8 

+2.9 

d 

+7.9 

-0.8 

+35.8 

+1,1 

e 

-2.9 

-9.2 

+24.7 

-5.0 

f 

-3.9 

-8.4 

+28 .5 

-5.0 

g 

-2.6 

-6.5 

+4.2 

-31-9 

h 

+1,7 

-3.6 

-7,5 

-21.2 

-10.3 

+1  .7 

—  i  r\  a 

IV  *  V 

-30.0 

* 

-Q  a 
♦  ** 

+4  ,G 

-12,4 

-29.3 

k 

-12.4 

->4.9 

-31,6 

+2 .1 

i 

-13.9 

+4.9 

—  ul  •  X 

+9.2 

m 

+0.5 

+  12.3 

-38.5 

+12.1 

n 

+1.4 

+13.0 

-44-6 

+12.2 

o 

+  3.2 

+14.8 

-44.3 

+42.0 

P 

+5.5 

+15. 3 

-35,1 

+32.0 

Sequence  of  sub -projectiles  listed  at  30  meters  is  not 
necessarily  identical  vcith  the  sequence  listed  at  10  meters , 
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Shot 

Impact  at 

10  Meters 

Impact  at 

30  Meters 

Number 

X(in.) 

Y (in. ) 

X{in.) 

Y (in, } 

7a 

+9.5 

-3.0 

+1.6 

+32 . 6 

b 

+6.0 

-3.0 

+6 . 5 

+33.8 

c 

+4.2 

-3.8 

+11.0 

+41  5 

d 

+0.9 

-10.3 

+29.5 

+1.9 

e 

-3 . 1 

-9,2 

+28.1 

-6.9 

F 

-4.0 

-9.3 

+17.0 

-9.5 

9 

-4.5 

-C.4 

+11.0 

-9.4 

h 

-9.7 

+1.0 

+24 , 6 

-34.7 

i 

-11.0 

+?.  8 

+  2 . 6 

-30.0 

j 

-10.9 

+3.9 

-15.8 

-24 ,0 

k 

-7.7 

+3.5 

-12.0 

-25,9 

1 

-1.3 

+9 . 8 

-25.8 

+  3.1 

m 

+1.4 

flO.l 

-33.8 

+8,4 

n 

+2 .6 

+12.3 

-36.  ) 

+10.4 

U 

+3.9 

-13.4 

-2  3,8 

+9.9 

P 

+10.1 

+0,5 

-4.0 

+30.2 

8a 

+7.0 

-4.7 

+23.3 

+20 . 3 

b 

+6,2 

-7.3 

+14.6 

+20,0 

c 

+3.7 

-7.0 

+16 . 0 

KL6.5 

G 

+1.0 

-7.4 

+15.6 

+13.0 

e 

-6.2 

-5 . 1 

+21.0 

-14.8 

f 

-7.8 

-3.2 

+19.  S 

-23.6 

Q 

-7.9 

-2,5 

+13.3 

-22.0 

h 

-7.9 

-2.3 

+3.4 

-25.5 

-4,6 

+8.9 

-19.6 

-1S.1 

3 

•5.0 

+9,2 

-23-4 

-14.6 

k 

-5.3 

+5.7 

-20.5 

-9,0 

X 

-5.9 

+10.2 

—211 

-2.0 

m 

+5,1 

+7.1 

-18.0 

+  30.6 

n 

+5  *  7 

+5.8 

-16,0 

+28.0 

o 

+5.6 

+5 . 5 

-13.7 

+24.3 

P 

+7,8 

+6 » 8 

-9.1 

+18.1 
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Shot 

Impact 

at  10  Meters 

Impact  at. 

30  Meters 

Humber 

X{in.) 

Y  {in. } 

X(in.) 

Y (in. ) 

9a 

+0.9 

+3.7 

+3.7 

+4.4 

b 

-1.9 

+1.9 

"7,1 

+5,1 

c 

-0.8 

+2.8 

+8.4 

+7.8 

& 

-4.6 

+2.8 

-2.8 

+8 ,6 

e 

-1.4 

+2.9 

-13.8 

+8.8 

T 

+2.2 

+3.0 

+7.9 

+  8.7 

g 

+0.4 

+3.1 

+8.2 

■r  9 » 0 

h 

+2  *  3 

+3.2 

+7.2 

+9.5 

i 

+2.0 

+3.3 

-3.7 

+9 . 7 

j 

+2.1 

+3 . 4 

+4  .5 

+10.3 

k 

+1.2 

+3.5 

+3.2 

+U. 6 

1 

-1.3 

+4.2 

+0.4 

+1.1 . 9 

Si 

-0.2 

+4.3 

-4.5 

+12.5 

n 

-2.1 

+5.2 

-4.3 

+15.7 

o 

+0.8 

+5.7 

+6.2 

X /  t7 

P 

+2.3 

+6.1 

+8.6 

i/A  4 
'  c  v 

10a 

-0.5 

-1.3 

+1.5 

-9,4 

b 

-0,4 

+0.3 

+2.8 

-3.7 

c 

0.0 

+1.9 

+2.8 

+1.5 

d 

-2.3 

+2.7 

-1.8 

+3.7 

e 

-1.0 

+2.8 

+7.3 

+4.4 

f 

+  0.8 

+2.9 

-3  7 

+4,6 

g 

-1.5 

+3.Q 

+6 . 9 

+4.8 

h 

-2,1 

+  3.1 

-3.4 

*4.9 

i 

-1.8 

+3.4 

-4.5 

+5.1 

-1-6 

+3 . 8 

+1.6 

+6,2 

k 

-1.3 

+3.5 

-l.S 

+6 ,5 

1 
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*  *  *  * 

+4.2 

-3.3 

+9,8 

m 

+4.8 

+2.7 
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n 

-3.5 

+5.6 

-4*7 

14.2 

o 

-2.2 
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«  i 
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V 

-1,4 
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+x5 . 3 

Sequence 

of  sub-pro je  utiles 

listed  at  30 

meters  is  a 

necessarily  identical  with  the  sc  :ence  listed  at  10  rasters 


APPENDIX  £ 


PAPEETE 3  ESTIMATION  BY  THE  METHOD  OF  MOMENTS 


A  a-1  -  h  -  f 

?(3)  x  1  X‘  * 


xi  for  each  i  -  X ,  x^iid 


The  expected  value  of  a  garuoa  distributed  random  variable 
is  a/A ,  the  variance  is  a/\2 .  To  use  the  method  of  moments 
see  the  sample  mean  equal  to  the  theoretical  mean  and  the 
sample  variance  equal  to  the  theoretical  variance.  Thus 

x  ~  a/X  and  5  2  -  c/X2 

Substituting  x  for  its  equivalent  expression  in  the  relation¬ 
ship  between  s:  and  the  distribution  parameters  vie id: 


sa  =  f  {*)  or  X  =  i- 

Substituting  the  value  for  lambda  in  the  expression  for  % 

—  n 

.  _  .  -  a  x  * 

y2.ei.es :  x  =  —  or  a  =  — r  - 

x  3 


Therefore,  the  method  of  moments  produces  the  following 
estimates  for  the  oaraaeters : 


Since  &  and  X  both  depend  On  x  and  s2,  they  ?rc  positively 
correlated. 

The  bias  b  (ai  of  an  estimate  is  defined  a£  b  to)  - 
E  (T)  -  9  where  T  rs  the  statistic  used  to  estimate  6,  the 
true  parameter  value.  Approximate  the  bias  by  using  a  first 
order  Taylor  expansion: 
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APPENDIX  C 

CENTER  OF  MASS  AS  A  FUNCTION  OF  RANGE 


The  information  about  the  center  of  mass  of  any  one  round 
is  limited.  The  location  of  ‘he  center  of  mass  is  known  at 
three  ranges:  0,  JLG,  and  3C  meters.  At  zero  range  the  center 
of  mass  is  at  the  origin  since  the  pellets  would  not  have  had 
an  opportunity  to  disperse.  The  location  at  10  and  30  meters 
can  be  calculated  from  the  data.  Nith  this  limited  informa¬ 
tion  it  is  difficult  to  accurately  depict  the  .cot ion  of  the 
center  of  mass .  Additional  experimentation  to  provide  more 
data  on  the  motion  of  the  center  of  mass  would  allow  a  sore 
ihcrcucrh  examination  of  the  notion  of  the  center  of  mass  as  a 


function  of  range.  A  desirable  situation  would  be  to  have 
sufficient  data  to  allow  an  analyst  to  determine  a  distribu¬ 
tion  of  the  center  of  mass.  In  the  absence  of  further 
information,  the  data  at  hand  should  be  used.  Therefore,  two 
possible  approaches  to  model  the  motion  of  the  center  of  mass 
are  presented  here.  Users  of  the  results  of  this  study 
should  determine  what  procedure  best  satisfies  their  view  of 
the  drift  of  the  center  of  mass. 

One  method  to  model  the  drift  of  the  center  of  mass  is 
for  each  one  of  the  coordinates  to  fit  the  lowest  degree 
polynomial  that  appears  reasonable.  An  examination  of  the 
data  would  seam  to  eliminate  a  constant  or  a  linear  relation¬ 


ship  between  either  coordinate  of  the  center  of  mass  and 
range.  Since  a  parabola  contains  three  parameters  and  coordinate 


values  are  available  at  three  ranees  for  each  shot  fire: 


it  is  possible  to  fit  a  parabola  through  the  coordinate 
values  for  each  coordinate  of  each  center  of  mass,  using  the 
parabolas  determined  for  the  two  coordinates  of  a  center  of 
sass,  the  position  of  the  center  of  mass  can  be  determined 
for  any  range.  However,  the  procedure  does  not  provide  any 


insight  into  the  physical  phenomenon  prt 
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the  notion. 


It  does  provide  a  method  of  modeling  the  motion  of  the 
center  of  mass  as  a  function  of  range  which  will  reproduce 
the  experimental  data. 

A  model  of  the  drift  of  the  center  of  mass  can  be  based 
on  the  physics  of  motion  if  one  can  accept  three  assumptions ; 
the  forces  acting  on  the  projectile  (the  recoil  and  air 
resistance)  act  fox  only  a  short  period  of  time,  once  the 
covering  on  the  pellets  is  aerodynamical iy  stripped  away 
air  resistance  is  negligible,  and  gravity  does  not  signi¬ 
ficantly  affect  the  pellets  c-vcr  their  maximum  effective 
range.  Under  these  assumptions,  the  forces  causing  the  center 
of  mass  to  drift  can  be  visualised  as  acting  instantaneously 
as  the  projectile  leaves  the  barrel.  After  leaving  the 
barrel,  the  projectile  would  travel  in  a  straight  line.  One 
can  then  use  the  two  data  points  at  10  and  30  meters  to 
estimate  the  parameter  of  a  straight  lina  through  the  origin 
using  the  least  squares  method.  The  line  mist  go  through 
the  origin  because  at  zero  range  the  center  of  mass  is  at 
the  origin. 


This,  z-odel  has  the-  appeal  of  being  based  on  physical 
principles  of  ration.  However,  it  has  the  disadvantage  of 
not  reproducing  the  actual  experimental  results  because  it 
Is  unlikely  that  a  least  squares  fitted  line  will  go  throat 
either  of  the  data  points. 

Once  the  user  has  selected  a  method  of  modeling  the 
action  of  the  center  of  sass  as  a  mr.cticr.  or  range,  ae  oa< 
use  that  model  to  generate  a  function  for  each  of  the  expo: 

rounds.  When  a  round  is  *  fired*  in  a  sinailation  of 


one  of  these  functions  based 


randomly  selected  to  provide  one  position  of  the  center  of 
isass  at  the  range  desired  for  the  simulation.  In  this  man 


the  drift  of  the  center  of  mss  can  be  reintroduced  to  the 


simulation. 


/  V  Vi 


ix,  ?}  ,  h  to,-  i2) 


-  S  CG.-V}  vsherti  the  variance  wztcx jL  V 


;i  agonal  and  has  the  f  ;m: 


VSEX  ^  VARX  -  2COV 


Inrs  the 


VhXZt-ViiSX- 2CGV  { Y ,  Y ) 
ataon  of  X  and  ¥  arc  identical ,  it 


possible  tc  determine  the  rntT3.es  n  the  V  jsafcri 
ting  the  COV  (x,Xi  teta. 


s  by  evaloa* 


COV  |X., xj  =  COV  fjt.iEX.) 

a  i 


Since  the  covariance  is  ?  linear  cps~$t<?i*t  the  constant  and 


the  Sira  car.  be  brought  outside  the  covariance  operator. 

CC¥  .  Aicsvc^.x.i  -  ii6..  -  I 


therefore,  the  V  matrix  has  descant t  l-~ 


Sc  tjr-X.Y-Y)  H{Q#{MrA-. 


Toe  square  or  too  distance  fr« 


ae  center  ot  sass  to  a 


pc-'- at  win;  aodt r led  coordinates 


(X-XI  tfy-v; 


{fill’! 

*  s  - 


Each  term  cn  the  right  side  of  the  last  equality  will  be 


recognized  as  the  square  of  a  standard  normal.  Hence 


After  1000  runs  through  the  program  ,  the  stored  test 
statistics  were  sorted  from  smallest  to  largest.  The  5,  10, 
25,  50,  75,  90,  and  95  sample  percentiles  were  printed.  The 
90  and  95  sample  percertiles  were  used  to  represent  the 
critical  values  for  this  goodness-of-fit  test.  The  sample 
percentiles  are  given  in  Table  3  as  a  function  of  the  alpha 
value . 


> 


TABLE  IV 

SAMPLE  PERCENTILES  FROM  GOODNESS -OF-FIT  SIMULATION 

SAMPLE  PERCENTILE 


ALPHA 


s 

10 

25 

50 

75 

on 

V  v/ 

95 

21.8 

24.5 

33,1 

84.5 

157.2 

1S5.3 

199. 

26,0 

31.2 

42.1 

65,2 

105.1 

135.5 

162. 

Cl 
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APPENDIX  F 


VARIATION  OF  LAMBDA  AS  A  FUNCTION  OF  RANGE 

Since  data  for  only  two  ranges  was  analyzed ,  only  two 
values  of  lambda  are  available  at  different  ranges.  An 
additional  value  for  lambda  can  be  deduced  by  observing  that 
the  variance  of  a  gamma  random  variable  is  a/A2  and  the 
variance  of  the  pellet  strikes  at  zero  range  is  essentially 
zero.  Since  a  is  a  constant  in  the  models  of  interest,  a 

must  increase  without  bound  as  the  range  decreases  to  zero. 

f 

i 

This  is  the  only  way  to  produce  a  zero  variance  at  zero 
range  within  this  family  of  distributions. 

Since  lambda  is  ur, bounded,  it  is  more  convenient  to 
determine  the  reciprocal  of  lambda  as  a  function  of  range 
and  require  that  function  to  go  through  the  origin.  This 
situation  bears  a  strong  reseaolence  to  the  problem  of 
determining  the  drift  of  the  center  of  r.ass  as  a  function 
of  range.  As  in  that  case,  it  would  hava  been  desirable  to 
have  additional  information  about  the  way  lambda  varied  as 
range  increased.  If  data  had  been  taken  every  five  or  ten 
meters,  a  more  detailed  analysis  of  the  variation  of  lambda 
could  have  been  accomplished. 

A  parabola  can  be  used  to  fit  the  three  values  of  the 
reciprocal  of  lambda.  Tnis  method  has  the  advantage  of 
exactly  duplicating  the  experimental  results  at  the  known 
ranges^  However,  it  offers  no  insight  into  the  reason  why 
lambda,  should  be  expected  to  vary  ;.n  this  way. 


An  alternative  method  can  be  cased  on  the  dimensionality 
of  the  reciprocal  of  lambda  and  the  moments  of  the  c-*mma 
distribution.  Since  lambda  varies  with  the  range  and  the 
reciprocal  of  lambda  times  a  constant  is  the  expected  value 
of  the  range/  then  the  reciprocal  of  lambda  must  have  dimen¬ 
sions  of  length.  Since  the  dimensions  of  the  reciprocal  of 
lambda  and  range  are  the  same,  this  suggests  that  a  linear 
relationship  between  the  reciprocal  of  lambda  and  range  might 
be  appropriate.  Therefore,  it  may  be  reasonable  to  use  the 
least  squares  technique  to  estimate  the  slope  of  a  line 
through  the  origin  and  use  that  line  to  determine  the  recipro 
cal  of  lambda,  and  hence  lambda,  as  a  function  of  the  range. 
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